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ABSTRACT 

The dust clouds in the torus of the quasar are irradiated by the central source, and the clouds at the inner 
radius of the torus re-radiate mostly in the near-infrared (NIR) wavebands. The ratio of the near-infrared 
luminosity to the bolometric luminosity LNiR/i'boi can therefore reflect the torus geometry to some extent. We 
find a significant correlation between the ratio of the near-infrared luminosity to the bolometric luminosity 
LNiR/i'boi and the central black hole mass Mbh for the Palomar-Green(PG) quasars, whereas no correlation is 
found between the Eddington ratio Lboi/i'Edd and the ratio LNiR/i'boi- Similar correlations are found for the 
mid-infrared and far-infrared cases. It may imply that the torus geometry, i.e., the solid angle subtended by the 
dust torus as seen from the central source, does not evolve with the accretion rate. The correlation of the solid 
angle subtended by the torus with the central black hole mass Mbh implies that the formation of the dust torus 
is likely regulated by the central black hole mass. We find that the torus thickness H increases with quasar 
bolometric luminosities, which is different with the constant torus thickness H with luminosity assumed in 
the receding torus model. The average relative thickness H/R of the tori in the PG quasars derived from the 
ratios of the infrared to bolometric luminosities is ~ 0.9. The further IR observations on a larger quasar sample 
including more fainter quasars by the Spitzer Space Telescope will help understand the physics of the dust tori 
in quasars. 

Subject headings: galaxies: active — quasars: general — accretion, accretion disks — ^black hole physics 



1. INTRODUCTION 

The putative dust torus is an important ingredient of the uni- 
fication model for active galactic nuclei (AGNs) (Antonucci 
119931) . A type 1 AGN spectrum is yielded if the broad-line 
region (BLR) is directly viewed through the hole of the face- 
on dust torus, while the BLR is obscured by the torus seen 
edge-on that leads to a type 2 AGN spectrum. This model 
can successfully unify Seyfert 1 and Seyfert 2 galaxies; ra- 
dio quasars and radio galaxies, etc (e.g., Antonucci & Miller, 
1985; Miller & Goodrich 1990; Tran 1995; Barthel 1989), 
which is supported by a variety of observational features of 
AGNs (e.g., Lawrence 1991; Hill, Goodrich, & Depoy, 1996). 

The dust torus is irradiated by the central engine of the 
AGN, and the dust cannot survive inside a critical radius 
at which the temperature is so high that the dust begins to 
sublimate. The temperature of the dust clouds at the in- 
ner edge of the torus is close to the sublimation tempera- 
ture. The irradiated dust produces thermal emission mostly 
at the NIR wavelength. The decomposition of spectral en- 
ergy distributions of AGNs suggested that the NIR con- 
tinuum emission (> 2/im) from Seyfert galaxies is domi- 
nated by thermal radiation from the hot dust surrounding the 
central engine (e.g., Kobayashi et al. 1993). The analy- 
sis on the infrared emission data of 64 Palomar-Green(PG) 
quasars favors the du st bein g heated by the central engines 
of AGNs (Haas et alJ l2003h . The inner radius of the dust 
torus is assumed to be roughly at the dust evaporation ra- 
dius j?(r,v.p) - 0.06(Lhni/10^^W )'/^ pc (Netzer & Laot 1996; 
iHill. Goodrich. & Depov' 1996). The intense monitoring ob- 
servations were carried out on the Seyfert 1 galaxy NGC4151, 
and a lag time of 48^3 days between the V and K light curves 
was measured. The inner radius of the dust torus in NGC415 1 
is ^ 0.04 pc corresponding to this measured time lag. For this 
source, the predicted evaporation radius /?(revap) — 0.015 pc 



(Netzer & Laor, 1996; Lboi ^ 6.5 x 10^*^ W, given by Kaspi et 
al., 2000), which is roughly consistent with the inner radius 
of the dust torus ^ 0.04 pc measured b y thermal dust rever- 
beration method llMinezaki et alJl2b04^ . If the radio quasars 
and galaxies are intrinsically same and can be unified by 
their different orientations, the relative thickness of the torus 
H /R ^ 2-3 is estimated from the fraction of quasars in the 
3CR sample (Hill. Goodrich. & Depoy 1996). 

If the relative thickness H/R of the torus is around unity, 
the random velocities should be >100km s"' for the matter 
in the tori at sub-parsec scales, because the dust torus is re- 
quired to be in dynamical equilibrium. If this motion is ther- 
mal, the corresponding temperature (^10^ K) is too high for 
dust to survive. It is therefore suggested that the cold dust is 
in the clouds moving at the random velocities > 100km s"' 
( Krolik & Be gelmanI Il988h . The origi n of such dust tori 
in AG Ns is still a mystery, though Zier & BiermannI (1200 IL 
l2002h suggested that the merger of binary black holes may 
lead the distribution of the surrounding matters to be a torus- 
like structure. In this paper, we perform statistic analysis on 
the properties of the dust tori in PG quasars. The cosmologi- 
cal parameters ^Im = 0.3, f^A = 0.7, and Hq = 70 km s~' Mpc"' 
have been adopted in this work. 

2. BLACK HOLE MASS 

'KasEi^LalJ ^od) derived an empirical relation between 
the BLR size and optical continuum luminosity for a sample 
of Seyfert 1 galaxies and quasars using the cosmological pa- 
rameters Hq = 75 km s"' Mpc"' and qo = 0.5, in which the 
sizes of BLRs are meas ured with the reverberation mapping 
method (lPetersonll993h . The relation 

ALa(5100)^"''"'±'"° 



/;blr = (26.4 ±3.0) 



h-days (1) 



lO'^'^ergs s" 

is derived bv iMcLure & JarvisI ( l200l for the same sample. 
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but the cosmological parameters used are same as this work: 
flu = 0.3, = 0.7, and Ho = 70 km s"' Mpc"'. 

The central black hole masses Mbh can be estimated from 
the velocities vblr of the clouds in the BLRs 



bh ■ 



^BLR^BLR 



(2) 



where the motions of the clouds are assumed to be virilized 
and isotropic. The velocities of the clouds in BLRs vblr are 
derived from the width of the broad emission lines. For most 
quasars, the BLR sizes have not been measured by the re- 
verberation mapping method, and the empirical relation Q is 
instead used to estimate the BLR sizes. The central black hole 
masses of quasars have been estimated from their broad-line 
widths and optical continuum luminosity (e.g., Laor, 2000; 
McLure & Dunlop 2001; Cao & Jiang 2002). For normal 
bright quasars, the bolometric luminosity can b e estimated 
from their optical luminosity La opt at 5100 by (iKaspi et al.l 
I2OOO,) 

Lboi — 9ALA,opt. (3) 

Assuming a constant accretion efficiency for all quasars, we 
have the conventional defined dimensionless accretion rate 



^Edd 



(4) 



3. TORUS GEOMETRY 



The dust clouds at the inner edge of the torus irradiated by 
the AGN central engine produce thermal emission mostly in 
the NIR waveband (e.g., Kobayashi et al., 1993), so the NIR 
luminosity of the dust torus is estimated by 



^^torus(^in) 

47r 



/inLbol ; 



(5) 



where Lboi is the bolometric luminosity of the AGN, Ai7toms 
is the solid angle subtended by the dust torus as seen from 
the central source at the inner radius Rin of the torus, and the 
factor /in describes the covering factor of the clouds at the 
inner radius of the torus and the uncertainties. The uncertain- 
ties arising from the NIR emission contributed by the central 
continuum emission of the AGN and starbursts in the AGN 
host galaxy may not be large, because this emission can be 
neglected compared with the NIR emission from the tori in 
most PG quasars ( IHaas et al. 2003). If the contribution in the 
NIR waveband from stars in the host galaxy can be neglected 
compared with that from the torus, the factor /„ describes the 
covering factor. Thus, we can infer the dust torus geometry 
from the NIR and bolometric luminosities of AGNs, 



Aaoms(«in) = 4^/i; 



-1 Lnir 

' Lbo\ ' 



(6) 



if the covering factor of the dust clouds at the inner edge of the 
torus is known. If we assume all the IR emission of quasars is 
from the dust torus irradiated by the AGNs, we can estimate 
the solid angle subtended by the whole dust torus by 



A(y^ =477 

'torus ^" 



Lho\ ' 



(7) 



where the covering factor / = 1 is adopted, because the puta- 
tive torus is required to be able to obscure the nuclear emis- 
sion while it is seen edge-on cAntonucci 199 3). If the torus is 
subtended at a constant soUd angle with radius (i.e., cone-like 



structure), the covering factor /„ at the inner radius .Rin can be 
estimated by 

/.„-^. (8) 

MR 

The solid angle estimated from Eq. (0 is an upper limit for 
some sources of which the IR emission from the host galaxies 
is important compared with that from the dust tori. 
The soUd angle subtended by the torus Ailtoms is 

AttH/R 



(4 + //2/7;2)l/2' 



(9) 



where H is the thickness of the torus. For the case of H/R< I, 
Equation can be approximated as 



An,, 



I-kH 
R 



(10) 



Combining the relations (|6j and (|9}, we can roughly estimate 
the torus thickness at its inner edge by 



, , _ 2/i„ j;in(LNIR/Lbol) 

"*^"'"^~rT — 77 — \2TT72 P ' 

l-/in (^NIR/Lbol)2 



(11) 



where the inner radius Ri^ of the torus is given by 
j^in - 0.06(Lhni/10^^W)'/2 pc (iNetzer & Laod 1199^ 
IHill. Goodrich. & DeDovll996l) . 



4. SAMPLE 



iHaas et aT] i2003l) provided a sample of 64 PG quasars with 
infrared spectral energy distributions (3- 150/im) observed 
by the ISO. The PG quasars have high infrared detection rate 
of more than 80 per cent. As discussed in Sect. 2, the cen- 
tral black hole masses in quasars can be estimated from their 
broad-line widths of H/3 and optical continuum luminosities 
ALa(5100). The full widths at half maximum (FWHM) of 
broad-line H/3 for 51 sources in this sample are available in 
iBoroson & Green ( 1992). For the remainder, they are at rela- 
tively high redshifts (z > 1). We search th e literature a nd find 
FWHM(H/3)=5100 km s"' for 0044+030 JBrothertodri996) ; 
and FWHM(H/3)=9601 km f\ 4613 km s"' f or 1634+706 
and 1718+481, respectively (Shields et al."2003"). It leads to a 
sample of 54 quasars with estimated black hole masses, which 
includes 41 radio-quiet quasars and 13 radio-loud quasars (4 
flat-spectrum and 9 steep-spectrum radio quasars). There are 
ten sources with FWHM(H/?)< 2000 km s"^ in this sam- 
ple (hereafter we refer to those quasars with FWHM(H/3)> 
2000 km s~' as broad-line(BL) quasars). 

5. RESULTS 

We estimate the central black hole masses of these PG 
quasars using their broad-line widths of H/9 and optical con- 
tinuum luminosities as described in Sect. 2. The NIR lumi- 
nosities Lntr of the P G quasars in the waveband of 3 - 10 //m 
given by iHaas et alJ 12003 ) are used in this work to explore 
the properties of the tori at their inner radii, because the emis- 
sion from the inner region of the dust torus irradiated by the 
AGN is dominant in the NIR waveband. This is supported 
by the t hermal dust reverberation measurements on NGC415 1 
iMinezaki et al. 2004) . In Fig. ^ we plot the relation between 
the ratios Lboi/^Edd and LNiR/Lboi- No correlation is found be- 
tween these two ratios. 

Haas et al. (2003) suggested a scheme of quasar evolution 
for their dust distribution surrounding the AGNs. The sources 
are divided into different classes according to their evolution- 
ary sequence. The classes and 1 represent cool and warm 
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ultra-luminous infrared galaxies respectively. Class 2 sources 
are young quasars, and the starbursts are still important in IR 
wavebands compared with the torus emission. There are two 
class 2 quasars in this sample, and we leave out these two 
infrared luminous quasars 0157 + 001 and 1351+640 in all 
our statistic analysis, because the contribution in the infrared 
waveband by the dust heated by the stars in these two quasars 
may not be neglected, though these two sources would only 
affect little on the statistic results. The relation between the 
black hole mass Mbh and the ratio LNiR/i^boi is plotted in 
Fig. 12] The generalized Kendall's r test (Astronomical Sur- 
vival Analysis, ASURV, Feigelson & Nelson 1985) shows 
that the correlation is significant at the 99.93 per cent level 
for the whole sample of 52 sources (the significant level is 
almost same at the 99.94 per cent for the sample of 42 BL 
quasars). The linear regression by parametric expectation- 
maximization (EM) Algorithm (ASURV) gives 

logioLNiR/Lboi = -0.097(±0.034)logioMbh/MQ-0.125(±0.279), 

(12) 

for the whole sample (the solid line in Fig. |2j, and 

logioLNiR/Lboi = -0.149(±0.038)logioMbh/MQ+0.325(±0.322), 

(13) 

for the broad-line quasars (the dotted line in Fig. |2j, respec- 
tively. Compared with the linear regression for the whole 
sample, the slope becomes slightly steeper for the BL quasars, 
while the correlation levels are almost same. For those 13 i 
radio-loud quasars in this sample, no statistic behavior sys- 
tematically different from the radio-quiet quasars is found. 

The BLR size Rblr is derived from the optical contin- 
uum luminosity at 5100 , so the estimated black hole mass 



covering factor /in at Rin can be estimated by using the rela- 
tion (|8}. Subtracting the three sources with LiR/Lboi ^ 1 and 
the infrared luminous source 1351 +640, we find that the ra- 
tios LNik/i'iR for the sources in our present sample are in the 
range 0.20-0.56 with an average of 0.39. We use Eq. il It 
to estimate the torus thickness at its inner radius Ri„ adopt- 
ing a covering factor /in = 0.39. The relations of Mbij-H(Rin) 
and L\,o\- H{Rin) are plotted in Figs. and |8] respectively. 
It is found that the torus thickness H at its inner radius Ri^ 
increases with the black hole mass Mbh and the bolometric 
luminosity Lboi. 

The solid angle Ai^toms subtended by the whole torus can 
be estimated by using Eq. Q from the ratio LiR/Lboi- We can 
then calculate the relative torus thickness H/R using Eq. (|9j. 
In Fig. |9] we plot the relation between the black hole mass 
Mbh and the relative torus thickness H/R. The three sources 
with LiR/Lboi Jc^ 1 do not appear in this figure. The statistic 
results are listed in Table 1 . 







-0.5 



Mbh oc VgL]^L°-^'(5100) is given by substituting Eq. Q into 
Eq. (|2}. The bolometric luminosity Lboi is estimated from the 
optical continuum luminosity by using the empirical relation 
(|3}. It is therefore doubtful that the correlation between the 
black hole mass Mbh and the ratio LNiR/i^boi may be caused 
by the common dependence of optical continuum luminosity 
Lx. We plot the relation between the bolometric luminosity 
Lboi and the ratio LNiR/Lboi in Fig. |3] It is found that the bolo- 
metric luminosity Lboi is correlated with the ratio LNiR/Lboi 
for the whole sample. In order to test the correlation between 
the black hole mass Mbh and the ratio LNiR/Lboi is affected to 
what extent by the common dependence of the optical contin- 
uum luminosity, we perform statistic analysis on the sources 
with 10"^^ erg s"' < Lboi < lO'*^'^ ergs"' (35 sources situated 
between the two dotted lines in Fig. |3}- We find that the ratio 
LNiR/Lboi is still significantly correlated with the black hole 
mass Mbh at the 97.25 per cent level for this sub-sample of 
35 quasars, while no correlation is found between Lboi and 
LNiR/Lboi (at the 75.01 per cent level). This means that the 
ratio LNiR/Lboi is intrinsically correlated with the black hole 
mass Mbh. 

We also perform statistic analysis on the relations: Mbh - 
LMiR/Lboi, Mbh-LpiR/Lboi, and Mbh- LiR /Lboi (see Figs, g] 
|6j. The mid-infrared (MIR, 10-40 /xm), far-infrared (FIR, 
40-150 /im), and total infrared (IR, 3 - 150 nm) luminosities 
of these quasars are taken from Haas et al. (2003). It is found 
that they are significantly correlated at the similar levels as the 
NIR case, while the slope becomes steeper for the FIR case. 
There are three sources in our sample have LiR/Lboi ^ 1 (see 
Fig. |6j. The detailed statistic results are summarized in Table 
1. 

The dust clouds at the inner radius Rin of the torus irra- 
diated by the AGN emit mostly in the NIR waveband. The 
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Fig. 1. — The relation between the ratios Lbol/iEdd and ^NlR/ibol- The 
circles represent radio-quiet quasars, while the squares and triangles represent 
steep-spectrum and flat-spectrum radio quasars, respectively. The sources 
with large circles have FWHM(H/3)<2000 km s"' . The relation between the 
ratios Lbol/i-Edd and LNiR/Lbol- 
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Fig. 2.— The relation between the black hole mass Mbh and the FiG. 4.— Same as Fig.lU but for the ratio Lmir/Li 

ratio Z-NlR/ibol- The solid line represents the linear regression for the 
whole sample but excluding two infrared luminous sources 0157 + 001 and 
1351 +640. The dotted line represents the linear regression for the sources 
with FWHM(H/3)> 2000 km s"' excluding two sources 0157 + 001 and 
1351 + 640. 
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Fig. 5. — Same as Fig.|2] but for the ratio LpiR/ibol- 



Fig. 3. — The relation between the ratios Lf,„i and LniR/^'bol- The dotted 
lines represent Lboi = 10*' erg s"' and 10** erg s"'. 
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Fig. 6.— Same as Fig.lH but for the ratio Lir(3 - 150Aim)/Lbol- ^^'O- The relation between the bolometric luminosity Lboi and the 

torus thickness H(Ria) at its inner radius i?i„. The solid line represents the 
linear regression for the whole sample but excluding the infrared luminous 
sources 0157+001 and 1351 +640. The dotted line represents the linear 
regression for the sources with FWHM(H/3)> 2000 km s"' excluding the 
sources 0157 + 001 and 1351+640. 
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Fig. 7. — The relation between the black hole mass M(,h and the torus 
thickness //(i?i„) at its inner radius The solid line represents the linear 
regression for the whole sample but excluding the infrared luminous sources 
0157 + 001 and 1351 +640. The dotted line represents the linear regression 
for the sources with FWHM(H/3)> 2000 km s"' excluding the sources 0157+ 
001 and 1351+640. 
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Fig. 9. — The relation between black hole mass Mbh and the relative thick- 
ness of the torus H/R. 
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6. DISCUSSION 

The NIR luminosities Ln tr of the PG quasars in the wave- 
band of 3 - 10 /im given by Haas et al. ( 2003) are used in this 
work. The NIR emission in the waveband > 2 /im is domi- 
nated by the thermal emission from the dust torus, which is 
supported by the decomposition of the SEDs of the Seyfert 
galaxies and the thermal dust reverberation measureme nts 
on NGC4151 (Koba vashiet al . 1993; Minezaki etliL 2004). 
The NIR emission ^ 2/i m is not inc luded in the NIR lumi- 
nosities Lnir provided bv iHaas et al ](|2003), so the solid an- 
gle f^torus subtended by the torus may be slightly larger than 
that estimated from the ratio LiMiR/i'boi by using the relation 
(|6j, No correlation is found between the ratios Lboi/^-Edd and 
LNiR/Lboi, which implies that the solid angle Ailtoius sub- 
tended by the dust torus does not evolv e with the accretion 
rate m. In the clumpy dust cloud model ( iKrolik & Begelmaiil 
the dimensionless torus thickness is estimated by 



where Vc is the typical random velocity of the dust clouds and 
Vk is the Keplerian velocity at radius R. This means the ra- 
tio of the typical dust cloud random velocity to the Keplerian 
velocity Vc/Vk is independent of the dimensionless accretion 
rate m. However, it should be cautious that the PG quasars are 
optically selected and they are very bright (most sources have 
Lboi/Lgdij > 0.1, see Fig. 0. We cannot rule out the presence 
of a correlation between the ratios Lboi/i'Edd and L^in/L\,o\, 
if more fainter quasars with lower values of Lboi/i^Edd are in- 
cluded in the statistic analysis. 

The intrinsic correlation between black hole mass Mbh and 
the ratio LNiR/iboi indicates that the dust torus geometry is re- 
lated with the central black hole mass. The larger black hole 
leads to a dust torus with smaller solid angle subtended by it 
seen from the central source. The black hole grows slowly 
during its bright quasar phase (usually less than a factor of 
2, see Kauffmann & Haehnelt, 2000), because the typical 
timescale for the bright quasar phase is ^ 10^"^ years (e.g., 
Yu & Tremaine, 2002). This means that the black hole mass 
estimated from the broad-line width and the optical contin- 
uum luminosity can roughly reflect the hole mass while the 
bright quasar was formed at an early time. If the dust torus 
was formed at nearly the same time as the bright quasar, it 
may imply that the formation of the dust torus is regulated 
by the central black hole, though the detailed mechanism is 
still unclear Zier & Biermann (2001, 2002) suggested a torus 
formation scenario that the merger of binary black holes may 
lead the distribution of the surrounding matter to be a torus- 
like structure. In this scenario, the black holes play important 
roles in the formation of the torus, which is in general consis- 
tent with our statistic results. 

Similar correlations between the black hole mass Mbh and 
the ratio LiR/Lboi are also found for the MIR and FIR cases 
(see Fig. 0]and|5}- This implies that the emission in the 
MIR and FIR wavebands is dominated by the thermal dust 
tori heated by the central AGNs, excepting the infrared lumi- 
nous quasars 0157 + 001. This source has very high ratios of 
the infrared to bolometric luminosities, especially in the MIR 
and FIR wavebands, which implies that the starbursts in this 
source contribute much in the MIR and FIR wavebands than 
the NIR waveband. It is still unclear why the slope of the 
correlation between Mbh and LpiR/i'boi is steeper than that for 
the NIR and MIR cases (see Table 1 for comparison). One 
likely explanation is that the torus properties (e.g., number of 



the clouds and their distribution, and the cloud optical depth, 
etc.) or/and the AGN spectral shape vary systematically with 
the black hole mass Mbh, because the spectrum emitted from 
the dust torus is determined by the torus properties, and the 
SED of the AGN (e.g., Nenkova, Ivezic, & Elitzur, 2002). 
The deta iled model calcul ations, for example, with the code 
DUSTY (iNenkova. Ivezic. & Elitzur 2002), are necessary for 
solving this issue, which are beyond the scope of this work. 
An alternative explanation is that the starbursts contribution in 
the FIR waveband is more important than the NIR and MIR 
wavebands, and the thermal emission from the dust torus irra- 
diated by the AGN contributes less in the FIR waveband. In 
this case, the starbursts contribution is supposed to change the 
slope of the correlation between Mbh and LpiR/i-boi, while it 
should still not be dominant to avoid destroying this correla- 
tion. This can be explained by recent IR observations showing 
that the nuclear starbursts are in the dusty tori around the nu- 
clei of some Seyfert galaxies flmanishi & Wadal2004ft . 

For most quasars in our sample, the IR emission is domi- 
nantly from the tori irradiated by the AGNs. In this case, the 
IR luminosity should be less than the bolometric luminosity. 
In Fig.|6l we find three sources with Lir(3- 150/im)/Lboi > 1. 
The IR emission from these four sources (including one class 
2 quasars: 0157+001; two narrow-line quasars: 0050+124 
and 1543+489) cannot be solely produced by the dust torus 
heated by the AGNs. The contribution of stars in the host 
galaxies should be important in these four quasars. The ratio 
LiR/Lboi ^ 1 can therefore be used as a criterion for starbursts 
dominant quasars. The detailed geometry of the dust torus is 
still unclear The opening angle of the torus may vary with 
radius. However, the ratio LiR/Lboi can always be used to 
estimate the solid angle subtended by the whole torus. For 
those four quasars with LiR/Lboi > 1, the IR emission is dom- 
inated by the starbursts and we are not able to estimate the 
solid angle subtended by the torus from the ratio LiR/Lboi- 
We find that the relative torus thickness H /R estimated from 
the ratio LiR/Lboi is in the range of 0.2-2.4 with an average 
of ~ 0.9 (see Fig. |9}. This is different from H/R^ 2-3 
estimated from the fraction of quasars in the 3CR sample 
jHill. G oodrich. & DeDov'1996V 

We may estimate the solid angle Ailtoms subtended by the 
torus at its inner radius from the ratio of the NIR to bolometric 
luminosities, because the NIR emission is dominated by the 
emission from the dust clouds at its inner radius irradiated by 
the AGN. The emission in the NIR waveband is less affected 
by the starbursts compared with the MIR and FIR wavebands. 
If the opening angle is constant for the torus at different ra- 
dius, the covering factor of the clouds at the inner radius can 
be estimated from the ratio L^i^/Lir (Eq. |8}. It implies that 
only a fraction of the nuclear radiation is absorbed by the dust 
clouds at the inner radius of the torus. The remainder may 
be absorbed by the dust clouds at larger radii, and they have 
lower temperatures than that at the inner radius. Most of their 
emission may be in the MIR or FIR wavebands depending on 
their distance from the nucleus. This is consistent with the fact 
that similar correlations are found between Mbh and LiR/Lboi 
for the MIR and FIR cases. We use a constant covering fac- 
tor fin = 0.39 to estimate the torus thickness H{Rin) at its in- 
ner radius. The slope of the correlation between Mbh-H(Rin) 
is ^ 0.23(see Table 1), which is obviously non-linear. The 
correlation between Lboi and H(Rin) has a slope of ~ 0.37, 
which is different from the constant torus thickness H with 
luminosity assumed in the reced ing torus model suggested by 
iHill. Goodrich. & DepovlJl996l) . The further IR observations 
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TABLE 1 
Results of statistic analysis 



Relation 


Sample 


Significant level (per cent) 


a 


h 


^>^bh ^NIR/^bol 


All 


99 93 


—0 nQ7 -1- 014 

\J .rjy 1 I 1)11 


_n 1 o-; -1- n 77Q 






99.94 


_() ]AQ + (1 038 


n ^25 + 322 


^bh ~ ^MIR / ^bol 


All 


99.94 


-0.1 35 ±0.040 


0.224 ±0.335 




BL quasars 


99.85 


-0.169 ±0.045 


0.526±0.381 


-^Wbh-ipIR/ibol 


All 


99.99 


-0.239 ±0.058 


0.658 ±0.475 




BL quasars 


99.88 


-0.255 ±0.073 


0.809 ±0.625 


Mbh-ilR/ibol 


All 


99.98 


-0.151 ±0.040 


0.697 ±0.330 


BL quasars 


99.93 


-0.178 ±0.044 


0.941 ±0.376 


Mbh-//m 


All 


99.51 


0.232 ±0.054 


-2.825 ±0.452 




BL quasars 


99.51 


0.275 ±0.061 


-3.220 ±0.523 


Lho\-Hm 


All 


> 99.99 


0.367 ±0.051 


-17.760 ±2.320 




BL quasars 


> 99.99 


0.344 ±0.044 


-16.708 ±2.019 




All(50) 


99.93 


-0.152 ±0.045 


I.Oll ±0.376 




BL quasars(42) 


99.93 


-0.210 ±0.054 


1.532 ±0.461 



Note. — Columns (3) and (4): coefficients of tlie linear regression: log^^^Y = alog^^^X +h. 



on a larger quasar sample including more fainter quasars by 
the Spitzer Space Telescope will help understand the physics 
of the dust tori in quasars. 

I am grateful to X.Y. Xia for stimulating discussion, and 
C. Zier for his helpful explanation of their torus formation 
model. 1 thank Caina Hao for pointing me an error in the 
original manuscript. This work is supported by the Na- 
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NKBRSF (grant G1999075403). This research has made use 
of the NASA/IPAC Extragalactic Database (NED), which is 
operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautic 
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